Spatial ordering of charge and spin in quasi one-dimensional Wigner molecules 
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Few-electron systems confined in quasi one-dimensional quantum dots are studied by the config- 
uration interaction approach. We consider the parity symmetry of states forming Wigner molecules 
in large quantum dots and find that for the spin-polarized Wigner molecules it strictly depends on 
the number of electrons. We investigate the spatial spin-ordering in the inner coordinates of the 
quantum system and conclude that for small dots it has a short-range character and results mainly 
from the Pauli exclusion principle while the Wigner crystallization in large dots is accompanied by 
spin ordering over the entire length of the dot. 



I. INTRODUCTION 

Strong confinement of charge carriers in two directions 
results in reduction of their degrees of freedom to a sin- 
gle one, i.e., in quasi one-dimensional motion. Such one- 
dimensional systems are realized typically in split-gate^ 
and cleaved-edge overgrowth— semiconductor quantum 
wires, as well as in carbon nanotubes^ but can also be 
realized in finite-size systems, i.e., in anisotropic quan- 
tum dots 5 or quantum rings^ There is a renewed in- 
terest in the one-dimensional systems related to the re- 
cent progress of vapour-liquid-solid fabrication of quan- 
tum wires of very high qualityiSiS 

The present paper is devoted to electron systems con- 
fined in one-dimensional quantum dots and in particu- 
lar to their Wigner crystallization 10 appearing when the 
electron-electron interaction dominates over the kinetic 
energy. Wigner electron solids (Wigner molecules) are 
predicted to appear in large dotsii or in strong mag- 
netic fields^ In the Wigner molecules the charge den- 
sity separates into distinct charge maxima each corre- 
sponding to one of confined electrons. Formation of 
Wigner molecules in the ground-state charge density in 
one-dimensional quantum dots was previously obtained 
in exact diagonalizationi^iiiii 5 . and density functional 
approach^ In one-dimensional dots the Wigner local- 
ization appears in the laboratory frame, in contrast to 
the inner-coordinate crystallization appearing in circu- 
lar quantum dotsji^ including quantum rings. Trans- 
port properties of Wigner crystals formed in open infi- 
nite one-dimensional systems have also been studiediiii^ 
The Luttinger liquid formalism has been applied^ to 
quantum wires with box- like boundary conditions, i.e., 
to the one-dimensional quantum dots. Melting of clas- 
sical one-dimensional Wigner crystals has recently been 
described^ 

We study the quasi one-dimensional quantum dots us- 
ing a configuration interaction approach with the effec- 
tive electron-electron interaction potential which we de- 
rived recently^! This work is a generalization of our ex- 
act two-electron studj^ to larger number of electrons. 
In the weak confinement limit the ground-state becomes 
nearly degenerate with respect to the spin configuration 



of the electron system^i 5 - Similar approximate degen- 
eracy has been found in quantum rings of large radius^ 
In this paper we study the parity symmetry of the nearly 
degenerate states forming Wigner molecules in large dots. 
We show that for spin-polarized electrons the Wigner lo- 
calization is formed only for one (even or odd) spatial 
parity of the state strictly dependent on the number of 
electrons. We present this dependence in the form of a 
theorem for which we provide a rigorous analytical proof. 
The found dependence of the parity of one-dimensional 
Wigner molecule states on the number of electrons is sim- 
ilar to the appearance of the magic angular momenta 
states for which Wigner crystallization is possible in cir- 
cular dotsi^iSi Furthermore, we discuss an inhibition of 
Wigner crystallization by a perturbation of the confine- 
ment potential through a central inversion-invariant po- 
tential well. 

Magnetic spin-ordering of electrons in one-dimensional 
space has been extensively studied 2 ^ in Hubbard mod- 
els which, in one dimension with only nearest-neighbor 
hopping interactions, predict the appearance of a 
low-spin ground stated This is a consequence^ of 
the Lieb-Mattis theorem 2 - 7 which implies that with- 
out spin-dependent interactions the ground-state of one- 
dimensional electron systems corresponds to the lowest 
possible spin quantum number (S=0 or 1/2). This fea- 
ture generally does not have to result in any spatial spin 
ordering. In this paper we use the exact numerical solu- 
tion of the Schrodinger equation to investigate the spatial 
distribution of spins in the one-dimensional quantum dot 
and the relation between the charge and spatial spin or- 
dering in the Wigner crystallization limit. We find that 
Wigner crystallization is accompanied by a long-range 
spin-ordering in the inner coordinates of the system in- 
stead of a spin-symmetry breaking predicted by density 
functional theorypiSiS In the ground-state this ordering 
has a clear antiferromagnetic character. 

This paper is organized as follows. In Section II we 
present the theoretical method. Section III contains the 
results for the Wigner localization and ground state de- 
generacy of the few electron systems. In Section IV we 
present proof for the dependence of the parity of spin- 
polarized Wigner molecules on the number of electrons. 
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Section V contains discussion of the effect of a central de- 
fect on Wigner crystallization. In Section VI the study 
of spin ordering is presented. Section VII contains our 
summary and conclusions. 



II. THEORY 

We consider N electrons confined in a quasi-one- 
dimensional quantum dot with strong lateral harmonic- 
oscillator confinement potential. The Hamiltonian of the 
system reads, 



N N N 



*=1 3>i 



where h stands for the single-electron Hamiltonian, 



2m* 



-V 2 



-(x 2 +y 2 ) + V(z), 



(1) 



(2) 



V{z) is the confinement potential in the z— direction. For 
a strong lateral harmonic-oscillator confinement energy 
(Hlu) the movement of electrons in the (x, y) plane is 
frozen to the harmonic-oscillator ground state. Then, 
one can perform integration^! over the lateral degrees of 
freedom which results in the following Hamiltonian, 

H = Nhu + Ef =1 h} D + Elt Ef> 4 (t/2)' («/0 

xerfc(z. iJ /2 1 /2/)exp(z2./2/ 2 ), (3) 

where Zij = \zi — Zj\ and 

h 2 d 2 



h 1D = 



2m* dz 2 



V(z) 



(4) 



is the single-electron one-dimensional Hamiltonian. In 
the following we will neglect the first term in Eq. 
i.e., the lateral confinement energy which is indepen- 
dent of the form of wave functions in the z direction. 
The last term in Eq. is the effective interaction 

energy2i for electrons in a quasi-one-dimensional envi- 
ronment resulting from integration of the Coulomb po- 
tential over the lateral coordinates, to* is the effective 
mass, k — e 2 /4ireoe, e is the dielectric constant, and 
I = yJhfmFuj. We assume V{z) — V we ii(z), a rectangu- 
lar potential well of depth Vb = 200 meV and width d. 
We adopt GaAs material parameters, i.e., to* = 0.067 
TO e o, e = 12.4 as well as fiu = 10 meV (I = 10.66 nm) for 
the lateral confinement energy. Calculations have been 
performed for TV = 2, . . . , 5 electrons by the configura- 
tion interaction approach with a basis set of Slater deter- 
minants built with single-electron spin-orbitals. Spatial 
single-electron wave functions have been obtained by nu- 
merical diagonalization of the finite-difference version of 
the single-electron one-dimensional Hamiltonian Q on 
a mesh of points. In construction of the Slater determi- 
nants with required spin and parity symmetries we use 



the spatial wave functions of up to 8 lowest-energy single- 
electron states which results in a Slater determinant basis 
size of up to 1520 elements and an accuracy better than 
0.01 meV. 

The present approach is based on the assumption that 
only the lowest state of the lateral (x, y) quantization is 
occupied. We performed test calculations for 2, 3 and 4 
electrons to check the validity of this approach. We al- 
lowed the electrons to occupy also the p-type lowest ex- 
cited state of the lateral quantization with angular mo- 
mentum ±h. Inclusion of p states not only allows for 
determination of the critical well length above which the 
p-shell is emptied, but it is also helpful to estimate the 
importance of the angular correlations in the x — y plane. 
The Coulomb matrix elements were evaluated using ef- 
fective interaction potentials derived with the use of the 
Fourier transform technique^! We have obtained the fol- 
lowing results: the p shell is left empty for d > 39 and 
41 nm for TV = 3 and 4, respectively (for 2 electrons 
the p shell is never occupied). Accounting for the x — y 
correlations via inclusion of the p-type orbitals in the 
configuration interaction basis lowers the 2-electron to- 
tal energy estimates by 0.18, 0.12, 0.01 and 10~ 4 meV 
for d = 40,50, 100 and 200 nm respectively. This "lat- 
eral correlation energy" value for the same values of d are 
equal to 0.4, 0.3, 0.08 and 4 x 10~ 3 meV for N = 3 and 
1.18, 0.67, 0.23 and 0.03 for N = 4, respectively. The 
energy overestimation in the range of d studied further is 
never significant and the present approach is nearly exact 
in the Wigner localization regime. 



III. GROUND STATE DEGENERACY AND 
WIGNER CRYSTALLIZATION 

In this paper we label the states by their total spin 
S and parity quantum numbers using the notation: S*±, 
where the positive (negative) sign stands for even (odd) 
parity. We discuss only the lowest-energy states for a 
given spin-orbital symmetry. Fig. ^a) shows the lowest 
energy levels of the 2-electron system multiplied by the 
dot length d as functions of d. For large dots the states 
0+ and 1— as well as 0— and 1+ become mutually de- 
generate. For large values of d potential energy related to 
penetration of electrons into the barrier region is negli- 
gible, the kinetic energy scales as 1/d 2 and the Coulomb 
energy as 1/d. Therefore, the product of energy and dot 
length for large d behaves as f(d) = C + D/d function, 
where the constants C and D are related to the Coulomb 
and kinetic energy, respectively. The energy levels of the 
degenerate pairs of states tend to different constants in 
the infinite d limit which is apparently due to different 
values of the Coulomb interaction in these pairs of states. 
The evolution of the charge density for growing length of 
the dot is shown in Figs. 1 (b-d). For large dots [cf. 
Fig. IHd)] the charge densities of the degenerate pairs of 
states become identical. In the ground state the charge 
density has two pronounced maxima which indicates the 
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d [nm] z [nm] z [nm] z [nm] 



FIG. 1: (a) Lowest energy levels multiplied by the dot length for N — 2. Numbers close to the curves denote the total spin 
quantum number of the corresponding states and signs +, — stand for even and odd parity symmetry, (b), (c), (d) - charge 
density of 0+, 1 — , 1+ and 0— states plotted with solid, dotted, dashed and dash-dotted lines for d — 50, 100 and 200 nm, 
respectively. 
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FIG. 2: (a) Lowest energy levels multiplied by the dot length for N = 3. (b), (c), (d) - charge density of 1/2—, 1/2+ and 3/2— 
states plotted with solid, dashed and dotted lines for d — 50, 150 and 200 nm, respectively. In (d) the charge density of the 
3/2+ state is shown by the dash-dotted curve. 



separation of electron charges into two charge islands, 
i.e., the Wigner crystallization. Fig. ^ shows that the 
singlet-triplet degeneracy obtained previously^ for the 
two-electron ground-state appears also in the first excited 
state. 

Fig. [2ta) shows the energy levels and Figs. HJb-d) 
the charge density for the lowest-energy states of the 
3-electron system for increasing d. For 3 electrons the 
Wigner molecule is formed in states 1/2—, 1/2+ and 
3/2— which become degenerate for large d. In the state 
3/2+ the charge density exhibits four maxima [cf. Fig. 
Efd)], which apparently prevents this state to be degen- 
erate with the ground state. 

In the 4-electron system the ground state corresponds 
to 0+ symmetry. The states 1 — , 1+ and 2+ for large 
dots [cf. Eta)] tend to the degeneracy with the ground 
state. The charge densities of these states for large dots 
present four distinct maxima [cf. Fig. Etd)]. Energy lev- 
els corresponding to states 0— and 2— are separated by 
a significant energy distance from the ground state [cf. 
Fig. 01a)] and in large dots they correspond to identi- 
cal charge densities with five maxima. The ground state 
charge density evolution obtained for N=3 and 4 is in a 



qualitative agreement with the results of Ref. £i£] . 

Finally, in the 5-electron system the ground- 
state of 1/2+ symmetry becomes degenerate with 
l/2-,3/2+,3/2-, and 5/2+ states [cf. Fig. Eta)] form- 
ing Wigner molecules for large dots [cf. Figs. Hfb-d)]. 
The spin polarized state of odd parity 5/2— does not 
become degenerate with the ground-state and its charge 
density in large dots forms six maxima [cf. Fig|3|(d)]. 

In the entire d range and for all electron numbers stud- 
ied the order of the lowest energy levels for given total 
spin quantum numbers (neglecting the parity) follow the 
order of the spin quantum numbers, which is in agree- 
ment with the theorem of Lieb and MattisiSi In large 
dots a ground-state degeneracy appears. In Ref. fii] 
the degeneracy was interpreted in terms of a vanishing 
tunnel coupling between the local minima of the total 
TV-dimensional potential energy. The present results in- 
dicate that the nearly degenerate states possess the same 
charge density in the laboratory frame. Moreover, we 
observe the following regularities. In the limit of Wigner 
localization the ground state of the TV-electron system 
appears for N different pairs of the spin and parity quan- 
tum numbers^* For even electron numbers N = 2 and 4, 
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FIG. 3: (a) 4-electron energy levels multiplied by the dot length, (b), (c), (d) - charge density of 0+, 1 — , 1+ and 2+ 4-electron 
states plotted with solid, dash-dotted, dotted and dashed lines for d = 100, 200 and 300 nm, respectively. In (d) the charge 
densities of 2— and 0+ states are marked with crosses and dots, respectively. 
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FIG. 4: (a) 5-electron energy levels multiplied by the dot length. Even (odd) parity levels are plotted with solid (dotted) 
lines, (b), (c) - charge density of 1/2+, 1/2—, 3/2+, 3/2— and 5/2+ states plotted with solid, dash-dotted, dotted, dashed 
and dash-double-dot lines for d = 100 and 200 nm, respectively. In (d) the charge density of the 1/2+, 3/2— and 5/2— state is 
shown by solid, dashed and dotted lines, respectively [charge densities of 1/2—, 3/2+ and 5/2+ are almost identical with the 
1/2+ and 3/2— charge densities are therefore omitted in (d) for the sake of clarity]. 



N charge maxima are formed only for even parity states 
with S = 0, while the odd parity zero-spin states possess 
N+l charge maxima [cf. Figs. 1(d) and 3(d)]. The spin- 
polarized Wigner-localized state can only be formed for 
one (even or odd) parity. Namely, the parity of the spin- 
polarized Wigner molecule state is even for 4 and 5 elec- 
trons and odd for 2 and 3 electrons. The charge density of 
the spin-polarized state of the other parity exhibits N + l 
maxima, i.e., the state does not form a Wigner molecule 
and as a consequence does not become degenerate with 
the ground state even for large dots. This conclusion will 
be cast into a theorem in the next section. 



IV. PARITY OF SPIN-POLARIZED WIGNER 
MOLECULE STATES 

Here we give an analytical proof of the theorem: for 
an odd number of electrons N — 2M + 1 as well as for 
an even number of electrons N = 2M the parity of one- 
dimensional spin-polarized Wigner-molecule state is even 
(odd) for even (odd) value of the integer M 



We will present here the proof for an odd number of 
electrons (the proof for even N can be easily deduced 
from the present demonstration) . For odd N one of the 
electrons resides near the center of the dot (point zq = 0, 
cf. Fig. |SJ|, and the others occupy spatially symmet- 
ric sites to the left and right of the dot around points 
which satisfy Zk = —z-k for k = ±1, ±2, . . . ± M. In the 
Wigner phase the total charge density possesses N max- 
ima corresponding to the separate single-electron charge 
densities. A single-electron charge | "lAfc ( ^) 1 2 is localized 
around point Zk- In the Wigner limit the overlap be- 
tween the single-electron charge densities vanishes (the 
proof is only valid when this overlap is negligible), so the 
total charge density can be expressed as their sum. Since 
the total charge density is symmetric with respect to the 
origin the following equality holds 

m-z)\ 2 = \^-k(z)\ 2 , (5) 

which results in the following relation for the single- 
electron wave functions 

M-z) = e iH i>-k{z), (6) 
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FIG. 5: Illustration to the proof that for N = 2M + 1 or 
iV = 2Af electrons the parity of the spin-polarized state which 
exhibits Wigner localization is accordant with the parity of 
M. 

where the phase 4>k is a real number. Relation © with 
changed sign of k reads 

i>- k {-z) = e*+-*Tl> h {z). (7) 

Phases 4>k and are not independent. Changing the 

I 



The phases in front of the determinant in Eq. 11411 cancel 
according to property (j5J). Exchanging M pairs of cor- 
responding columns in the determinant we arrive at Eq. 
(|12|l but multiplied by (— 1) , which proofs that the par- 
ity of spin-polarized one-dimensional Wigner molecule 
state is determined by the odd or even value of M. 

We have found that 2- and 4-electron zero-spin states 



sign of z in Eq. Q and making use of relation J6j we 
arrive at 

^ k (z) = e'*-*^-*) = e^+^ty-kC*), (8) 

hence 

4>k = -<t>-k, (9) 

up to an unimportant multiple of 2-7r. Considering rela- 
tion © for k — and reminding that we arrive at the 
same value V'o(O) (non-zero for odd N) approaching the 
origin from both positive and negative sides we arrive at 
<^o=0 and consequently tpo is an even parity function 

M-z) = M*)- ( 10 ) 

Since the considered state is spin-polarized the spin and 
spatial parts of the wave function can be separated into 
a product 

X{zi,<Ji, ■ ■ ■ , z N ,a N ) = a{ai)a(a 2 ) . . . a(a N ) x 

X$(zi,z 2 ,...,z N ), (11) 

where a is an eigenfunction of the single-electron spin 
z-component operator. The spatial wave function W can 
be written as a Slater determinants^ 



(12) 



(13) 



(14) 



can form a Wigner-localized charge density only for even 
spatial parity. We are unable to proof in general that the 
zero-spin state with Wigner localization has to be of even 
parity for even N. But for TV = 2 such a proof is easily 
given. In this case the spin and spatial parts of the wave 
function can be separated as follows 

X 0+ {z\,cri,Z2,cr2) = [a((Ti)/3((T 2 ) - a(a 2 )/3(<Ti)] x 



^(z 1 ,z 2 ,...,z N ) 



1p- M (zi) 1p-M+l(zi) 



.lpM-l(z 2 ) 4>m(z 2 ) 



iJJ-m(zn) 1p-M+l(z N ) ...lp M _i(z N ) 1pM(ZN) 

We apply the parity operator on \I/ and make use of properties (6) and (9) obtaining 



V(-Zu-Z2, 



-ZN 



e-^ M i> M {zi) e-^'-^M-i(zi) 
e-^ M ip M {z 2 ) e-t+v-^M-ifa) 



^-^-M+lizi) e^V- M (zi) 



e-^ M ip M (z N ) e-^'-^ M -i(z N ) ...e**"->-il>-M+i(z N ) e*"il>- M (z N ) 
Phase factors can be extracted from each of the determinant columns, which yields 



V(-Zl,-Z 2 , 



-zn) 



= g — i($M+<i>M-l + ---+<i>-M+l+<t>-M) 



^Pm(zi) ^m-i(zi) . . .ip-M+i(zi) V-m(zi) 
^m(z 2 ) ip M -i{z 2 ) . . -4>-m+i(z 2 ) 4>_ M (z 2 ) 



4>m{zn) iPm-i(zn) 

I 



.iP-m+i{zn) iP-m{zn) 
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x [ih(zi)il>-i(z2) + ip-i(zi)i>x(z 2 )] • (15) 

Applying the parity operator to the spatial part of this 
wave function and making use of the properties of the 
single-electron wave functions given above we find that 
this wave function is of even parity. Moreover, it fol- 
lows that construction of a symmetric spatial wave func- 
tion for odd-parity singlet (zero-spin) two-electron states 
(0— ) requires at least three single-electron functions, for 
instance, the function 

^°~(zi, z 2 ) = ipo(zi)ih(z2) + ^1(21)^0(22) 

—ipo(zi)ip-i(z2) - ^-1(^1)^0(^2), (16) 

is of odd parity provided that we take zero phase shifts 
in relation ©. Indeed, the 0— state for N = 2 exhibits 
three charge maxima [see Fig. 1(d)]. Moreover, con- 
struction of a triplet antisymmetric spatial wave func- 
tion with even parity (1+) also requires at least three 
localized functions, for instance 

^ 1+ (ZI,Z 2 ) = 1po(zi)lpl{z 2 ) - Tpl(zi)->p (z2) 

+ipo(zi)ijj-i(z 2 ) - ip-i(zi)ipo(Vi), (17) 

possess the required symmetries for zero phase shifts in 
relation JfJJl. The charge density corresponding to wave 
functions l|ltj|) and (|17|l is the same provided that the 
overlaps between the functions ipi are negligible. Fig. 
1(d) shows that the charge densities of the states 0— and 
1+ are indeed indistinguishable. The area below the cen- 
tral maximum of the probability density of degenerate 0- 
and 1+ states in Fig. 1(d) is two times larger than the 
area below each of the extreme maxima, which can be 
interpreted by saying that one of the electrons stays in 
the neighborhood of the center of the system with 100% 
probability while probabilities of finding the other one at 
the left or right end of the well are equal to 50%. This fea- 
ture is in agreement with the probability amplitudes 
and H7J. Although in the wave functions (|16fl and i|17|) 
the electron positions are separated, this separation has 
a non-classical character since the charge maxima at the 
left and right ends of the dot correspond to sub-electron 
charges. Therefore we do not refer to this separation as 
Wigner localization. Average electron-electron distances 
in states described by wave functions and H17fl are 
smaller than in states 0+, 1— with two charge maxima, 
which leads to a larger value of the Coulomb interaction 
energy and consequently to an energy separation between 
pairs of degenerate states 0+,l— and 0—, 1+ presented 
in Fig. ^a) in the weak confinement limit. 



V. WIGNER CRYSTALLIZATION IN THE 
PRESENCE OF A DEFECT POTENTIAL 

The presence of defects can significantly perturb the 
Wigner crystallization in large systems. We consider here 
a thin attractive cavity just deep enough to bind one 



electron. The perturbed quantum dot potential is of the 
form 

V(z) = V weU (z) + V defect {z), (18) 

where Vd e fect{z) = —50 meV for \z\ < 1 nm and 
Vdefect(z) = for \z\ > 1 nm. The assumption that the 
defect is localized in the center of the system does not 
perturb the inversion invariance of the total potential. 
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FIG. 6: (a) Lowest energy levels for N — 2 as functions of the 
length of the well with a central attractive cavity [Eq. USU I. 
Numbers close to the curves denote the total spin quantum 
number of the corresponding states and signs + , — stand for 
even and odd parity symmetry, respectively, (b), (c), (d) 
- charge density of 0+, 1 — , 1+ and 0— states plotted with 
solid, dotted, dashed and dash-dot curves for d — 50, 100 and 
200 nm, respectively. 

Fig. EJa) shows that contrary to the unperturbed 
quantum well potential [cf. Fig. QJa)] the 0— and 1+ 
states become degenerate with the 0+ and 1— states. 
Figs. HHb-d) show the evolution of the charge densities 
of the four considered states with increasing size of the 
system. For large well thickness [cf. Fig. [Jd)] the charge 
densities of these states become indistinguishable. One 
of the electrons is trapped by the potential of the cen- 
tral cavity which results in the sharp central peak of the 
charge density. The probability to find the other electron 
at the left or right side of the origin are equal. This differs 
essentially from the two-electron Wigner molecule charge 
density in the unperturbed dot [cf. Fig. ^ , for which the 
probability to find an electron in the center of the well 
was negligible and for which each of the two charge max- 
ima could be associated with an integer electron charge. 
The formation of three maxima in the charge density is 
possible for all states [cf. Eqs (|16|) and l|17fl . for 0— and 
1+ states, similar formulas can be given for the other 
two]. Therefore, the ground-state tends to a fourfold de- 
generacy in contrast to the double degeneracy for the 
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unperturbed dot [cf. Fig. 1(a) 
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not perturb the number of charge maxima, Wigner local- 
ization appears similarly as for the unperturbed dot [cf. 
Fig. |2] for 1/2+, 1/2—, and 3/2— states which become 
degenerate in the Wigner localization limit. State 3/2+, 
which according to the theorem given in Section IV can- 
not form a Wigner phase lies higher in energy, like for 
the unperturbed dot. 

The influence of the central attractive defect is quali- 
tatively different for odd and even electron number. For 
an odd number of electrons it simply enhances the lo- 
calization of the central electron, and does not influence 
the ground-state degeneracy. While for even N it de- 
stroys Wigner crystallization leading to the appearance of 
an extra charge maximum corresponding to sub-electron 
charge and allows more states to become degenerate with 
the ground-state. 



VI. SPATIAL SPIN ORDERING IN THE 
WIGNER LIMIT 



FIG. 7: (a) Lowest energy levels for N = 3 as function of the 
length of the well with a central attractive cavity [Eq. 1181 1. 
Numbers close to the curves denote the total spin quantum 
number of the corresponding states and signs +, — stand for 
even and odd parity symmetry. (b),(c),(d) - charge density of 
0+, 1 — , 1+ and 0— states plotted with solid, dotted, dashed 
and dash-dot curves for d = 50, 100 and 200 nm, respectively. 

Fig. shows the lowest energy levels and correspond- 
ing charge density evolution for the 3-electron system. 
Contrary to the 2-electron system the central defect does 



It is interesting to look whether the low-spin ground 
states exhibit any spatial antiferromagnetic ordering of 
the electron spins. For even number of electrons and ar- 
bitrary dot length the ground-state corresponds to zero 
total spin. In this case the spin-up and spin-down densi- 
ties are exactly equal to each other, so that spin-ordering 
is not visible in the laboratory frame of reference. In or- 
der to investigate a possible spin-ordering in the zero-spin 
ground states one has to look into the inner coordinates 
of the quantum system. We use here the spin-dependent 
pair correlation functions (PCFs) defined for a given state 
by the expectation values 



N N 



F s P a ^{z a ,z b ) = {Y.Y. 5 ^ - Zi ^ Zb - z i) (W ff X^) >< a (°i) a (°j)\ + \P(<n)P{<rj) >< P{^)P{a 3 )\) ) , (19) 
\i=i j>i I 



and 



N N 



F° p p P° F (z ai z b ) = (Y^T, 6 ^ - z ^ Zb - z j) (K^fo) >< a i°i)P{°j)\ + >< P(<Ti)a(tTj)\) j , (20) 

\i=l j>i I 



where a and (3 stand for spin-up and spin-down eigen- 
states respectively. Functions (I19I20|I give the probabil- 
ity of finding at positions z a and z b a pair of electrons 
with the same <|19[) or opposite (|20|l spins. The sum of 
functions (|19fl and H2U|) gives the spin-independent PCF. 

Fig. IHa) shows the PCF plots for the 4-electron 
ground state in a small quantum dot [cf. Fig. ffib)] 
with d — 100 nm. The position of one of the electrons is 
fixed near the right end of the dot [position marked by 



the thin vertical line in Fig. IHJa)]. We see that the prob- 
ability of finding an electron with the same spin in the 
neighborhood of the fixed-position electron is zero, which 
is a signature of the Pauli exclusion principle. At the left 
side of the dot probabilities of finding an electron with 
the same or opposite spin as the one of the fixed position 
electron are nearly equal. For the total zero-spin states 
in relatively small dots the spin-ordering in the inner co- 
ordinates is of short range and results from the Pauli 
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FIG. 8: Pair correlation functions (PCFs) for 4 electrons in state 0+ (a,b) and state 1— with S z = h (c) for d — 100 (a) and 300 
nm (b,c). One of the electrons is fixed and its position is marked by a thin vertical line. Solid curves show the spin-independent 
PCF, dashed (dotted) curves show the opposite (same) spin PCF. 



exclusion. We only found a long-range inner-coordinate 
spin-ordering in the Wigner crystallization limit. Fig. 
IHtb) shows the plot for the 4-electron ground state with 
d = 300 nm. The charge density of the system exhibits 
4 distinct maxima [cf. Fig. |3fd)]. We fix the position of 
one of the electrons at the rightmost density maximum 
[cf. the vertical line in Fig. Etb)]. The probability that 
the electron in the adjacent maximum has the opposite 
spin is nearly 100 %. The spin-dependent PCFs also dif- 
fer for the two charge maxima at the left of the origin. 
An electron confined at the first (second) charge maxi- 
mum to the left of the origin is more probable to have the 
same (opposite) spin as the one of the fixed electron. The 
ordering is of a probabilistic character, so that the anti- 
ferromagnetic order of spins is the most probable to be 
found, but the probability is not 100%. The spin-ordering 
in this state has a clearly antiferromagnetic character and 
its range covers the entire length of the dot. A similar 
inner-coordinate antiferromagnetic order was previously 
found for quantum rings. 6 

The 100% probability of finding the opposite spin in 
the charge maximum adjacent to the maximum associ- 
ated with the fixed electron presented in Fig. Efb) is not, 
as one could naively expect, related to the Pauli exclu- 
sion. In Fig. [Hfc) we plotted the PCF for the 1— state, 
which becomes degenerate with the ground 0+ state in 
the weak confinement limit. We see that in this state the 
spin of electrons confined in the two central maxima is 
independent of the spin of the electron at the rightmost 
maximum. However, in this state one may expect that 
the electrons at the opposite ends of the dot have the 
same spin, which means that also in this state a long- 
range spin-ordering exists, even if it is not of antiferro- 
magnetic origin. 

For odd number of electrons the difference between 
spin-up and spin-down densities appears in the labora- 
tory frame. This is qualitatively different from quantum 
rings, which in fact are endless structures. Fig. EJa) 
shows the spin densities for a relatively small dot length 
of d = 100 nm [too small for the ground-state Wigner 



localization to appear, cf. Fig. 0[b)]. The spin-up elec- 
trons tend to gather at the extreme left and right ends 
of the dot as well as in its center. The spin-down density 
is minimal in the center of the dot, and the overall spin 
density (difference of the spin-up and spin-down densi- 
ties) exhibits antiferromagnetic sign oscillations within 
the dot. These sign oscillations are due to the electron- 
electron interaction since in the noninteracting electron 
system the majority spin-up density is nowhere smaller 
than the spin-down density. For larger systems [d = 250 
nm, cf. Fig. EJb)] the antiferromagnetic spin oscillations 
become more pronounced. However, for even larger d [cf. 
Fig. EJc) and (d)] , for which the Wigner molecule appear 
in the 1/2+ ground state, the typically antiferromagnetic 
real-space spin-ordering with the spin orientation chang- 
ing between the adjacent charge maxima vanishes. 

Let us look at the spin distribution in the inner coor- 
dinates of the 5-electron 1/2+ ground state. Fig. HOf a) 
shows the PCF plots for d = 100 nm. Electrons of the 
same spin as the fixed electron do not appear in its close 
neighborhood, but are more probable to be found at the 
center of the dot than electrons of opposite spin. Prob- 
ability of finding an electron at the opposite side of the 
dot is independent of its spin. The spin order in this 
relatively small dot (d = 100 nm) is clearly short-range 
which is similar as for the case of 4 electrons in a small 
dot [cf. Fig. Ufa)]. The PCF plots for opposite spins at 
the left end of the dot start to differentiate for d = 200 
nm [cf. Fig. Efb)]. For d — 300 nm, for which Wigner lo- 
calization is observed [cf. Fig. 4(d)], the PCF plots show 
a long-range antiferromagnetic spin ordering. Notice the 
growth of the PCF plot for the same spin direction in the 
closest neighborhood of the fixed-position electron from 
d = 200 to 300 nm in Figs. Efb) and (c). Pauli exclusion 
plays a less significant role for larger distances between 
the charge maxima. 

Density-functional studiesi 6 ^^ predict the appearance 
of interlocked waves of opposite spins in the laboratory 
frame for long quasi one-dimensional dots. The appear- 
ance of the spin-density wave for even electron number 
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FIG. 9: Spin-up (solid lines) and spin-down (dashed lines) 
densities for the ground-state 5-electron system 1/2+ with 
S z =h/2 for different system sizes. 



amounts in spin symmetry breaking. Recently^ it was 
found that for even N the formation of the spin density 
wave in the density functional theory accompanies the 
Wigner crystallization. But in the present study we find 
that for the exact solution spin symmetry is conserved 
and Wigner crystallization is associated with the inner 
space spin ordering. In the exact solution the interlocked 
spin densities in the laboratory frame can only be ob- 
served for odd numbers of electrons, but the presented 
5-electron case shows that this effect is not necessarily 
related with Wigner crystallization. In the exact solu- 
tion the electrons with opposite spins avoid one another 
in the inner space. A mean field approach can only ac- 
count for this effect by symmetry breaking. The reason 
for the occurrence of spin symmetry breaking in the mean 
field approach for large single-dimensional dots are sim- 
ilar to the origin of the broken spatial symmetry mean 
field solutions for the magnetic field induced Wigner crys- 
tallization in circular structures^ 

In large systems the spin-independent PCF plots be- 
come identical for all states degenerate with the ground- 
state [cf. spin-independent PCFs for the 4-electron de- 
generate 0+ and 1— states in Figs. IBJb) and (c)]. This 
means, that in Wigner-molecule states electrons avoid 
one another with the same efficiency independently of 
their spins. As a matter of fact this is the origin of the 
appearance of the ground-state degeneracy in the Wigner 
molecule regime. One-dimensional Wigner molecules 
present pronounced magnetic properties related to the 
long-range spin-ordering in the inner coordinate space. 
This ordering for different degenerate spin eigenstates 



may be typical for ferromagnetic, antiferromagnetic or 
even an other form of order. Due to the vanishing energy 
spacing between the different spin states the spin mag- 
netic properties of Wigner molecules are of a very soft 
character. The Wigner molecules should be extremely 
susceptible to any spin-dependent interactions. In par- 
ticular even a weak additional effect promoting the spin- 
polarized phase can result in spin-polarization of the sys- 
tem. A possible spin-polarization of the one-dimensional 
electron gas has been found? in transport measurements. 



VII. CONCLUSIONS AND SUMMARY 

We have studied the ground and excited states of elec- 
tron systems confined in quasi-one-dimensional quantum 
dots using an exact diagonalization approach. For large 
systems we found Wigner localization which appears not 
only in the ground state but also for several excited states 
which eventually leads to the degeneracy of the ground- 
state in the large d limit. We have considered spin and 
spatial parity of states forming Wigner molecules. We 
have shown that the parity of the spin-polarized state 
which forms a Wigner molecule is strictly determined by 
the number of electrons. 

We have discussed the effect of a central attractive de- 
fect which destroys Wigner crystallization for an even 
number of electrons allowing more states to become de- 
generate with the ground-state in the weak confinement 
limit. For odd electron numbers the central defect en- 
hances the localization of the electron occupying the cen- 
tral position in the Wigner molecule and does not affect 
the ground-state degeneracy. 

We have investigated the spin ordering effects associ- 
ated with Wigner crystallization. We have found that for 
small dots the spatial spin-ordering in the inner coordi- 
nates has a short-range character and results mainly from 
the Pauli exclusion principle. The long-range spatial spin 
order appears only in the Wigner molecule regime when 
the electrons occupy distinct sites within the quantum 
dot. We conclude that in one-dimensional quantum dots 
the Wigner crystallization is a necessary condition for 
the long range spin ordering to appear. We have identi- 
fied the effect of spin symmetry breaking observed in the 
density functional theory as a tendency of the mean field 
method to mimic the internal-space spin ordering present 
in the exact solution for the Wigner molecule regime. 
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FIG. 10: Pair correlation functions (PCFs) for 5 electrons in the state 1/2+ for d = 100 (a), 200 nm (b) and 300 nm (c). The 
position of the fixed electron is marked by a thin vertical line. Solid lines show the spin-independent PCF, dashed (dotted) 
lines show the opposite (same) spin PCF. 
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The approximate ground-state degeneracy of a N electron 
system in large one-dimensional dot equals 2^ (taking into 
account the degeneracies with respect to the ^-component 
of the total spin) as reported already in Ref. £— ]. For 4 and 
5 electrons also higher-energy states for a given spin and 
parity quantum numbers (not discussed in the present pa- 
per) become degenerate with the ground state. We have de- 
termined the parity of these states. Namely, the 4-electron 
ground state becomes degenerate with second 1— and 0+ 
states and the 5-electron ground state becomes degener- 
ate with the second 1/2—, 3/2—, 1/2+ and 3/2+ symmetry 
states as well as with the third 1/2+ state. 
A single determinant is sufficient to describe the spin- 
polarized Wigner molecule formed by N electrons local- 
izated around N space sites. See also the discussion of the 
electron-electron correlation given in Ref. £—]. 



